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Hopf transformation, which enable us to analyze the asymptotic behaviors of general 
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Finally, we discuss the generalization to a one-dimensional conservative neighborhood- 
n CA, i.e., PCAn. 
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1. Introduction 

Cellular automata (CA) are dynamic systems in which space and time are discrete 
and physical quantities take on a finite set of discrete values. Despite their simple 
construction, CA exhibit complicated behavior and generate complex patterns. By 
virtue of these properties, CA have been used as mathematical models for complex 
phenomena in the fields of physics, chemistry, biology, economics, and sociology[Tl |2]. 

Recently, one- dimensional conservative CA have attracted much attention as 
models for traffic fiow[3]. One of the simplest models is Rule 184, which is a 
neighborhood-three CA (also called an elementary CA). Since Rule 184 exhibits phase 
transition from a free-fiow state to a congestion state as real traffic fiow does, it has 
been used as a basic model and extended to many other CA models. However, despite 
the intensive research on CA models for traffic flow, there are few studies dealing with 
all one-dimensional conservative CA in a unifled way. 

A powerful method for dealing with CA is the ultradiscretization methodjH [5], 
which connects difference equations and CA by applying the following simple formula 
for the transformation of variables: 



A conservative CA called the "box and ball system" (BBS) is obtained from the discrete 
KdV equation using this method, and other CA that possess a solitonic nature are 
obtained from soliton equations. Moreover, Rule 184 is derived from the discrete Burgers 
equation [6]. 

Equation ([1]) shows that the addition of the difference equations corresponds to the 
max operation of the ultradiscrete equation. The multiplication and division correspond 
to + and — , respectively, since we have 



Though the subtraction is not well deflned for the ultradiscrete equation, we can 
automatically obtain the ultradiscrete equation and its solution by replacing +, x and 
/ by max, + and — , respectively, if the subtraction is not included explicitly in the 
difference equations and their solutions. The algebra defined by operations of max, +, 
and — is called "max-plus" algebra. Binary operations such as AND, OR, and NOT 
can be easily converted to max-plus operations. 

In a previous study[7], we investigated one-dimensional neighborhood-four 
conservative CA, referred to as particle cellular automata four (PCA4), and found that 
evolution equations for PCA4-1, 4-2, and 4-3 can be obtained in the form of max-plus 
expressions in combination with "min," i.e., "max-min-plus" expressions, which are the 
key to solving the equations and analyzing asymptotic behaviors. 

In this study, we investigate PCA5s and show that the evolution equations for 
PCA5s that belong to certain types can be obtained in the form of max-min-plus 
expressions from fundamental diagrams. In addition, we show that the equations 




(1) 



£log(e^/" X e^/") =A + B, 



e\og (e^/7e^/^) =A~B 
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obtained are transformed into other max-min-plus expressions by ultradiscrete Cole- 
Hopf transformation, which enables us to analyze the asymptotic behaviors of general 
solutions. Furthermore, we give the Lagrange representation for the class of PCA5s, 
which describes the motion of particles and can be obtained from a fundamental diagram. 
Finally, we discuss the generalization to the neighborhood-n CA case, that is, PCAn. 



2. Particle CA 



Let us consider the following equation 



u 



n+l 

'j 



f{u]_i, u]_i^„ u]^,) i-oo < J < oo, < n), (2) 



where n and j are integers representing the timestep and space site number, respectively, 
and I and r are positive integer constants. Assuming that u takes the value of or 1, ([2]) 
is an evolution equation for a one-dimensional CA with the rule defined by R{= l + r + 1) 
variable function /. We call this CA a " neighborhood- i? CA." Following Wolf r am [H [2], 
to each Rule /, we assign a rule number N{f) such that 

{ui,u2,...,uit)e{o,i}^ 

In this study, we focus on the one- dimensional CAs that satisfy 

K K 

5Z ^ XI ■/'("i-^' '"i-^+l, ■ • • 5 Uj+r) (3) 

where K{> R) is the size of the one-dimensional lattices, and mj+a' = Uj holds for all j. 
From ([3]), the condition 

x:«-=i:<-^\ (4) 

i=i j=i 

holds, which means that the sum of u, that is the number of Is at all space sites, is 
conserved for arbitrary n. Let denote the number of particles at the jth site and 
the nth timestep, and each "1" in the solution represents a particle. Then, particles 
move among sites according to the evolution rule defined by (jj]) without creation or 
annihilation. We call the CA satisfying this condition (jl]) a "particle CA" (PCA) in 
this sense. 



3. PCA5 



In [8], Hattori and Takesue showed that a neighborhood- i? CA with Rule / is a PCA if 
and only if / satisfies 

U2, ur) - ui+i = q{ui, U2, M/j-i) - q{u2, Ms, ... , un), (5) 



q{ui, U2, Ur-i) = ^Mfe - ^/(p, 0, ■ • ■ ,0, Ui, U2, UR_k), (6) 
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for all U2, ■ ■ ■ , Ur) G {0, 1}-^. 

Applying Hattori and Takesue's result to neighborhood-five CA expressed by the 
equation 

we obtain a system of 2^ linear equations with 2^ variables /(a, 6, c, d, e) for all 
{a,b,c,d,e) G {0,1}^ from ([5]) and (E]). By finding binary solutions to the system, 
we obtain 428 rules for one-dimensional conservative neighborhood-five CA (PCA5). 
Reflection symmetry 

fi{a,b,c,d,e) = /2(e, ci, c, 6, a) 

and Boolean conjugation symmetry 

/i(a, b,c,d,e) = 1 — /2(1 — a, 1 — 6, 1 — c, 1 — rf, 1 — e) 

for rules /i and /2 reduce the 428 rules to 129 rules. Furthermore, by excluding the 
rules for neighborhood-three and -four CA (PCA3 and PCA4, respectively), we have 
the following 115 rules: 



2863377064 
2945232012 
3099375756 
3103688898 
3136326348 
3153627912 
3167259896 
3180117376 
3200290508 
3213933712 
3218115792 
3220734008 
3370699752 
3421555648 
3450663328 
3484741764 
3488923844 
3705199640 
3771264248 
4040228048 



2881005752 
2945363074 
3099506818 
3103754432 
3136457410 
3163077816 
3167521996 
3182211488 
3200487104 
3214980000 
3218639896 
3220996108 
3372009676 
3422079752 
3452751296 
3486043912 
3489185848 
3705721824 
3822120144 



2881267852 
2945428608 
3099572352 
3116153216 
3136522944 
3163470978 
3167653058 
3184299456 
3201011208 
3216027824 
3219162080 
3221127170 
3384605056 
3431529656 
3454845408 
3486567632 
3639663552 
3706245928 
3822644248 



2881398914 
2947326124, 
3102247912 
3120335296 
3137047048 
3163536512 
3167718592 
3186393568 
3202384620 
3216289924, 
3219686184 
3366517672 
3388787136 
3431791756 
3464560268 
3487091736 
3640187656 
3706769648 
3824214256, 



2881464448 
2947457186, 
3103295736 
3132144268 
3148921728 
3165434028 
3169616108 
3196108428 
3202581216, 
3217067968 
3220209904 
3367565496 
3400596108 
3435711736 
3482385552 
3487613920 
3703627712 
3707031748, 
3824738360 



2944969912, 
3098065832, 
3103557836, 
3132275330, 
3153103808, 
3165630624, 
3169812704, 
3198202540, 
3203105320, 
3217592072, 
3220472004, 
3367827596, 
3404778188, 
3448569216, 
3484479664, 
3488138024, 
3704151816, 
3707293752, 
3888178416, 



which is the smallest rule number among their equivalent rules. In this paper, we assign 
a number "m" in the range 1-115 to each rule, and call them PCA5-m. 
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4. Obtaining evolution equations in the form of max-min-plus expressions 
from a fundamental diagram 

In the field of traffic-fiow theory, a fundamental diagram is a useful tool for determining 
the traffic state of a roadway. This diagram relates traffic ffux and traffic density. Under 
periodic boundary conditions, the density is defined by 

1 ^ 

i=i 

where K is the period of sites. Since the number of particles is conserved, p is a constant 
irrespective of time n. The average fiux is defined by 

For large enough n, the evolution often approaches a steady state where (p becomes a 
constant. In case of convergence, the constant is defined by 

Q = lim g". 

n— >oo 

The graph of Q versus p is called a "fundamental diagram." 

The following subsection shows that the evolution equations for PCA5s that belong 
to certain types can be obtained in the form of max-min-plus expressions from a 
fundamental diagram. 

4.1. Type- A 

Let us take PCA5-34 (Rule 3163536512) as an example. The evolution equation for 
PCA5-34 is given by 



together with the following rule table of /. 



abode 


11111 


11110 


11101 


11100 


11011 


11010 


11001 


11000 


10111 


10110 


10101 


10100 


10011 


10010 


10001 


10000 


f{a, b, c, d, e) 


1 





1 


1 


1 


1 








1 











1 


1 


1 


1 


abode 


01111 


OHIO 


01101 


01100 


01011 


01010 


01001 


01000 


00111 


00110 


00101 


00100 


00011 


00010 


00001 


00000 


f{a, b, c, d, e) 


1 





1 


1 


1 


1 








1 
























The upper row of the table shows all possible combinations of binary arguments, and 
the lower row gives the value of /. 

From ([5]) and ([6]), the evolution equation for PCA5-34 can be rewritten as 

together with the following rule table of q. 



abed 


1111 


1110 


1101 


1100 


1011 


1010 


1001 


1000 


0111 


0110 


0101 


0100 


0011 


0010 


0001 


0000 


q{a, b, c, d) 





1 


1 


1 








1 


1 





1 


1 


1 
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1/31/2 2/3 1 
P 




p = 0.2 p = 0.4 p = 0.6 p = 0.8 

Figure 1. Fundamental diagram and space-time evolutions of PCA5-34 



The upper row of the table shows all possible combinations of binary arguments, and 
the lower row gives the value of q. 

Employing numerical simulation for ([8]), we obtain the fundamental diagram and 
examples of evolutions for PCA5-34 shown in Figure [1] 

Since the fundamental diagram in Figure [1] is a piecewise linear curve, the 
relationship between Q and p can be expressed in terms of the max-min-plus expression, 

Q{p) = max(min(2p, 1 — p), min(p, 2 — 2p)), (9) 

which is a piecewise linear function composed of four hnear functions 2p, 1 — p, p, and 
2-2p. (See Figured) 

Let us consider replacing variables mp in 



by 



mp 



Then, we obtain 



k=l 



k=l 



(m > 0) 
+k-i < 0) 



(10) 



Q(m^_2, u], u]+^) = max(min(u^_2 + u]_^, 1 - u]),mm(u]_^, 2 - u] - u]^^)) (11) 
from Q. (See Figure [31 ) 
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Figure 2. Four linear functions of p and the composed piecewise linear function 




Figure 3. Obtaining expression for flux from the fundamental diagram 



Since Q{u^_2, u^-i, uf, "u^+i) takes exactly same values as those of the rule table for 
g{u'j_2, u^_i, u'j , ujj^i) of dH]) for all possible combinations of binary values, (|H]) can be 
rewritten as 

= u] + max(min(M"_2 + u]_^, l-u]), mm{u]^^, 2 - u] - u]^^)) 

- max(min(M"_i + m", 1 - min(M^, 2 - - M^+g))- (12) 

This is the max-min-plus expression for PCA5-34. Introducing the ultradiscrete Cole- 
Hopf transformation from u to F, 

= Ff- F^_„ (13) 

we obtain the evolution equation in the form of max-min-plus expressions for F, 

= min(max(F,"_2, i^^+i - 1), max(F,-_i, - 2)), (14) 

which is composed of linear functions of -FJYfc (^ ^ ^)- Being able to express evolution 
equations for F as a combination of linear functions of F is critical for analyzing the 
asymptotic behavior of general solutions of PCA. We used this fact for the analysis of 
the asymptotic behavior of PCA4 in a previous study [7]. 
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Note here that fll4p can be directly obtained from (jH]) by applying the following 
replacements. (See Figure HI) 

mp + a — )■ -F"_^ — a 
max — )■ min 
min — )■ max 




Figure 4. Obtaining equations for F from the fundamental diagram 



It is known that PCAs allow two different representations: i^w/er representation and 
Lagrange representation. In the Euler representation, particles are observed at a certain 
fixed point in space as dependent (field) variables, while in the Lagrange representation, 
we trace each particle and follow its trajectory. Thus, a dependent variable represents 
each particle's position in the Lagrange representation. 

In previous studies[9l [10], we proposed an Euler-Lagrange transformation for 
PCAs by developing the following transformation formulas for the variable of Euler 
representation m", which denotes the number of particles at the jth site and nth 
timestep, and the variable of Lagrange representation xf, which denotes the position of 
the ith particle at the nth timestep, 

u] = F^- Ff_,, (15) 

N 
i=l 

where H{x) is the step function defined by H{x) = 1 if x > 0, and H{x) = otherwise. 
Applying the transformation f|T6|l to f|T4l) and using the formulas, 

N N N 

HU - min(ai, k)) = max(J2 ^(j - «^), 5Z ' ^*))' ^^^^ 

i=l k=l k=l 

N N N 

J2 H{] - max(a,, 6,)) = min(5^ H{j - a,), H{j - 6,)), (18) 

i=l k=l k=l 
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max(^ H{j - ai) - m, 0) = ^ H{j - ai+„), (19) 

i i 

where we assume that ai < a2 < ■ ■ ■ < aN and 6i < 62 < • ■ ■ < ^at, we obtain 

x1+' = max(min(x^ + 2, x^+i - 1), min(x^ + 1, - 2)), (20) 

which is the Lagrange representation for PCA5-34 in the form of a max-min-plus 
expression. Note here that ( 120|) can be directly obtained from (|9]) by applying the 
following replacements. (See Figure |5l) 

mp + a^Xi+a + ^ (21) 




1 
Figure 5. Obtaining equations for x from the fundamental diagram 



To summarize, the procedure we have performed above is as follows. 

(i) Employ numerical simulation for PCA and obtain the fundamental diagram. 

(ii) If the fundamental diagram is a piecewise linear curve, construct Q{p) from the 
fundamental diagram. 

(iii) Construct flux q from Q{p) and obtain the evolution equation (Euler 
representation), 

(iv) Obtain the evolution equation for F, 

= ^iF^_„ F-_,, F^, i^'Vi, i^%). 

(v) Obtain the evolution equation for x (Lagrange representation), 

, Xj — ij Xj, Xj-\-i, Xjjf.2)- 

Among the 115 rules for PCA5, there are 17 for which the evolution equations 
for u", F", and x" can be obtained in the form of max-min-plus expressions from 
fundamental diagrams by using the procedure described above. These equations are 
given in Appendix A. In the tables of the Appendix, m denotes the number of PCA5-m 
and N{f) denotes the rule number defined by Wolfram. Hereafter, we call the 17 rules 
type- A. 
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Fundamental diagram 

piecewise linear function Q{p) 

\. Replacing variables from p to u 

max-min-plus expression for flux 
q{uj-2,Uj-i,Uj,Uj+i) 

Ultradiscrete 
I Cole-Hopf 
^ Transformation 

max-min-plus expression for F 
4>{Fj-2, -Fj-i, Fj, Fj+i, Fj+2) 

\^ Euler-Lagrange transformation 
max-min-plus expression for x 

Figure 6. Procedure for obtaining evolution equations of type-A PCA5 in max-min- 
plus expression 



12. Type-B 

Other than type-A, there exist PCA5 rules for which the fundamental diagram is a 
piecewise linear curve, and the evolution equations for m", F", and can be obtained. 

Let us consider PCA5-15 (rule 3099572352). Employing numerical simulation, we 
obtain the fundamental diagram for PCA5-15, as shown in figure [3 figure [7] indicates 
that the piecewise linear curve can be composed from two linear functions 2p and 2 — 2p 
as follows. 

g(p) =min(2p, 2- 2p) (22) 
Following the same procedure as in the case of type-A, we obtain a flux as follows: 

q{uj^2, Uj-i, Uj, Uj+i) = min(Mj_2 + 2 - Uj - Wj+i) (23) 



Evaluating the values of fl2^ for all possible combinations of binary values, we obtain 
the following table. 



abed 


1111 


1110 


1101 


1100 


1011 


1010 


1001 


1000 


0111 


0110 


0101 


0100 


0011 


0010 


0001 


0000 


q{a,b, c, d) 





1 


1 


2 





1 


1 


1 





1 


1 


1 














However, the flux of PCA5-15 is g 


iven 


by 




















abed 


1111 


1110 


1101 


1100 


1011 


1010 


1001 


1000 


0111 


0110 


0101 


0100 


0011 


0010 


0001 


0000 


q{a,b,e, d) 





1 


1 


1 





1 


1 


1 





1 


1 


1 
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and the value of g(l, 1,0,0) is found to be different. Thus, (1231) is not the flux of 
PCA5-15. 

To solve this problem, let us suppose that the piecewise linear curve in Figure [7] is 
not composed of two but more than two linear functions. Let us consider 

g(p) =min(l,2p,2-2p) (24) 

instead of fl22|) . Note here that fj22l) and fl2^ are identical functions for real values of p. 
(See Figure El) 

2i 1 




1/2 
Figure 7. Fundamental diagram of PCA5-15 





Figure 8. mm{2p, 2 - 2p) and min(l, 2p, 2 - 2p) 



If we apply variable replacement from p to m to the above expression, we obtain 
q{uj-2, Uj_i, Uj, Uj+i) = min(l, Uj_2 + Uj-i, 2 - - Wj+i) (25) 
which gives the following table: 



abed 


1111 


1110 


1101 


1100 


1011 


1010 


1001 


1000 


0111 


0110 


0101 


0100 


0011 


0010 


0001 


0000 


q{a,b,c, d) 





1 


1 


1 





1 


1 


1 





1 


1 


1 
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which is identical to the flux of PCA5-15. Thus, the evolution equation in the form of 
a max-min-plus expression for PCA5-15 is 

= + min(l, u]_, + u]_,, 2 - u] - u]^,) - min(l, u]_, + u], 2 - u]^, - u]_,,).{26) 
We can also obtain the evolution equations for F and x as follows. 

= max(i^" - 1, - 2) (27) 

^ = min(xr+i, + 2, x^^, -2)=x^ + min(2, x^Vi - x^, x^^^ - a:? - 2) (28) 

There are nine rules for PCA5, of which the evolution equations for u, F, and x are 
obtained using the procedure described above. These equations are given in Appendix 
B. We call the nine rules type-B. 

4-3. PCA5 rules other than type- A and type-B 

We have employed numerical simulations for all 115 PCA5 rules to obtain fundamental 
diagrams, from which we have obtained evolution equations in the form of max-min-plus 
expressions for 17 rules of type-A and nine rules of type-B. There are 89 rules for PCA5, 
for whom evolution equations have not yet been obtained. From the results of numerical 
simulations, the 89 rules are classified into the following two cases. 

• Although the obtained fundamental diagram is a piecewise linear curve, replacing 
variables from p to u does not give the correct flux expression for q. In addition, 
the type-B procedure does not seem to work well. 

• As the obtained fundamental diagram is not a piecewise linear curve, we cannot 
use the procedure of replacing variables from p to u. 

At this point, it is not clear if there may be other procedures for obtaining max-min-plus 
expressions for them. 

5. PCAn 

We can extend our analysis in the previous section to neighborhood-n CA, i.e., the 
PC An case. 

Let us consider the case of n = 6. One example of PCA6 is Rule 
13755053124876288240. From numerical simulation, we obtain the fundamental diagram 
for the rule shown in Figure [H Following the procedure of type-A of PCA5, we obtain 

q{uj-2, Uj-i, Uj, Mj+i, Uj+2) = max(— Uj, min(Mj__2 + Mj-i — 1,1 — Uj — Uj+i — Uj+2)), 

min(-Uj_i — 1,1 — Uj — Mj+i, min(Mj_2 + — 2,2 — Uj — Uj+i — Uj+2)), 

(29) 

= mm{F^^„ max(i^"_2 + 1, F]"^, - 1), max(i^"_i + 1, - 1), max(i^'!_2 + 2, F]"^, 
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-1 







1 



Figure 9. Fundamental diagram for Rule 13755053124876288240 



x^+i = max(xl^ - l,mm(x^i + 2,a;I!,i - 3),min(x^i + l,x^+i - 2),mm(x^2 + 2,x^+2 - 



Although we have not investigated all PCA6 rules, the analysis of PCA5 in the previous 
section and the example of PCA6 above assure us that our procedure can be applied to 
PCAn rules. 

6. Summary 

We have studied PCA5 and shown that the evolution equations for type-A and type- 
B can be obtained in the form of max-min-plus expressions from a fundamental 
diagram. The obtained equations have been transformed into max-min-plus expressions 
that are composed of linear functions of F by ultradiscrete Cole-Hopf transformation. 
Furthermore, we have obtained the Lagrange representation of the evolution equations. 

Although we have not obtained evolution equations for all 115 PCA5 rules, it is 
important for us to have been able to introduce a unified approach using max-plus 
algebra, together with the fundamental diagram, for examining PCAs. 

In a previous study [7], we analyzed asymptotic behaviors of solutions and derived 
functions Q{p) for PCA4 mathematically. We have not done this for PCAS, but will be 
able to do so starting from the equations for F obtained in this study. This problem 
will be addressed in a future study. 

Finally, investigating the generalization to the neighborhood-n case is important 
for the future, and we will report on this in a forthcoming paper. 
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Appendix A. Type-A 



m 




q{uj-2,Uj-i,Uj,Uj+i) 


i 


zodoo77(Jd4 


min {uj-2 + Uj-i, 1 — Uj — Uj+i) 


o 
O 


0001 OCJTOCO 


max {—Uj, min(?ij_2 + Uj-i — 1, 1 — uj — uj+i)) 


4 




max {—Uj — Mj+i, min(uj_2 + Uj-i — 1,1 — Uj — Uj+i)) 


r 




OOO 1 A£i A A AO 


min (wj-2 + Uj-i, 2 — Uj — Wj+i) 


QQ 
OO 




max {—Uj — Uj+i, min (uj-2 + Uj-i — 1, —Uj) , min {uj—i — 1, 1 — Uj — Uj+i)) 


O A 


oiooooooiz 


max (min (uj-2 + Uj-i, 1 — Uj) , min {uj-i, 2 — Uj — Uj+i)) 


QC 
OO 


olo /ozlyyo 


max {—Uj, min (?ij-2 + Wj-i — 2, 1 — m^- — Uj-^-i)) 


on 


31d7d5305o 


max {—Uj — Uj+i, mm {Uj-2 + Uj-i — 1, —Uj) , mm {Uj-2 + Uj-i — 2,1 — Uj — Wj+i)) 


4U 


3167 1 i859z 


max (min {uj-2 + uj-i, 1 — uj) , min (w,j-2 + uj-i — 1,2 — Uj — Uj-^.l)) 


50 


3200487104 


max (min (uj-2 + Uj-i, 1 — Uj — Uj+i) , min (uj-i, 1 — Uj) , min (wj-2 + Uj-i — 1,2 — Uj — Mj+i)) 


54 


3203105320 


max (min {uj-2 + uj-i, 1 — uj — Uj+i) , min (uj-2 + w^-i — 1, 2 — m^- — m^+i)) 


59 


3217067968 


max (min (wj-i, 1 — uj) , min {uj-2 + uj-i — 1,2 — uj — 


64 


3219686184 


max (min (mj-i , 1 — Uj — itj+i) , min (Mj-2 + Wj-i — 1, 2 — itj — Uj+i)) 


65 


3220209904 


max (0, min (uj-2 + Wj-i — 1, 2 — — Wj+i)) 


68 


3220996108 


max {—Uj, min {uj-2 + — 1: 2 — — 


69 


3221127170 


max [ — 11 J — iij+L-j mill ("j-2 + '^j-i ^ 1- 2 — f j — "j+i jj 


96 


3488138024 


max(min(uj_i, 1 — — Uj+i), min(uj_2 + Uj-i — 1, 1 — Uj)) 






TO 


Nif) 


(j){Fj-2 , Fj-i, Fj , Fj+i , Fj+2) 




1 


2863377064 


max(Fj_2,J>+2 - 1) 


3 


2881267852 


min(Fj+i,max(Fj_2 + l,Fj+2 - 1)) 


4 


2881398914 


min (F,+2, max {Fj-2 + 1, Fj+2 - 1)) 


5 


2881464448 


max(F,_2.-F,+2 -2) 


33 


3163470978 


miu (_rj-+2. max {Fj^2 + 1. ^y+i) • m-i-^ + 1. - Ijj 


34 


3163536512 


min (max (J}_2, - 1) , max {Fj-i,Fj+2 — 2)) 


38 


3167521996 


min(F,-+i,max(F,_2 + 2,Fj+2 - 1)) 


39 


3167653058 


min {Fj+2, max (Fj_2 + 1, Fj+i) , max {Fj-2 + 2, Fj+2 - 1)) 


40 


3167718592 


min (max (-Fj-2, Fj+i - 1) , max (-Fj-2 + 1, Fj+2 - 2)) 


50 


3200487104 


min (max (-Fj_2 , -Fj+2 - 1) , max (Fj-i , Fj+i - 1) , max (Fj-2 + 1, -Fj+2 — 2)) 


54 


3203105320 


min (max (Fj -2, Fj+2 - 1) ,max(Fj-2 + l:Fj+2 - 2)) 


59 


3217067968 


min (max (Fj-i, Fj+i - 1) ,max(Fj-2 + l,Fj+2 — 2)) 


64 


3219686184 


min (max (Fj-i, Fj+2 - 1) , max (Fj-2 + 1, Fj+2 - 2)) 


65 


3220209904 


min (Fj, max (Fj_2 + 1, Fj+2 - 2)) 


68 


3220996108 


min (Fj+i, max (F,_2 + 1- Fj+2 - 2)) 


69 


3221127170 


min (Fj+2, max (Fj_2 + 1, Fj+2 - 2)) 


96 


3488138024 


min (max (Fj-i, Fj+2 - 1) , max (Fj-2 + 1, Fj+i - 1)) 
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m 








i 




min(xi + 2, x^+i — 2) 


Q 

o 


OQQ1 O^J'TQKO 


max(a;i — 1, min(xi_i + 2, Xi+i — 2)) 


4 


OQQ1 Q0Q01 A 


max(a;j — 2, min(a;,_i + 2, Xi+\ — 2)) 





zooi40444o 


min(.Tj + 2, — 2) 


QQ 
OO 




max(a;i — 2, min(xi_i + 2,Xi — 1), min(a;i_i + 1, ccj+i — 2)) 


o4 


C)iDc50c5DOiZ 


max(min(a;i + 2, a;j+i — 1), min(xi + 1, Xi+2 — 2)) 


OO 


oiD / oziyyo 


iTiax(.T^ — l,niin(xi— 2 ~t~ 2^2',i^i — 2)) 


oy 


'Jit; Tt; n o 
olD / DOoUoo 


max(.Tj - 2, min(xi_i + 2, x^ - 1), min(a;j_2 + 2, Xi+i — 2)) 


/I n 


olo i ( looyz 


max(min(a;j + 2, — 1), niin(a;i_i + 2, a;i+2 — 2)) 


OU 


oZUU4o / iU4 


niax;(min(a;j + 2, — 2), min(a;j + 1, ajj+i — 1), min(a;i_i + 2, a;j+2 — 2)) 


04 


ozUoiUoozU 


max(min(xi + 2, — 2), min(a;i_i + 2, Xi-1-2 — 2)) 


59 


3217067968 


max(min(xi + 1, Xi+i — 1), miii(a;i_i + 2, Xi+2 — 2)) 


64 


3219686184 


max(min(a;j + 1, Xj+i — 2), min(a;i_i + 2, Xi-^-2 — 2)) 


65 


3220209904 


niax(0, min(xj_i + 2, a;j+2 — 2)) 


68 


3220996108 


max(xi — l,inin(a;j_i + 2,Xi+2 - 2)) 


69 


3221127170 


max(xi — 2,min(x,;_i + 2,.t,;+2 — 2)) 


96 


3488138024 


max(min(xi + 1, x^+i — 2), min(a;i_i + 2, x^+i — 1)) 


Appendix B. Type-B 




m 


N{f) 


g(%-2,%-l,Mj,Mj + l) 


2 


2881005752 


mill (max (0, Mj-2 + — 1) , 1 — — Uj+i) 


13 


3099375756 


max {—Uj, min (0, Mj-2 + Uj-i — 1, 1 — Uj — 


14 


3099506818 


max {—Uj — Wj+i, min (0, Uj-2 + Uj-i — 1,1 — Uj — Uj+i)) 


15 


3099572352 


min (1, Uj-2 + Uj-i,2 — Uj — Uj+\) 


53 


3202581216 


max (min {uj-2 + Uj-i, max (0, 1 — Uj — Uj+i)) , min {uj-2 + Uj-i — 1, 2 — — Wj+i)) 


63 


3219162080 


max (min (uj-i, max (0, 1 — uj — Uj^i)) , min (uj_2 + w^-i — 1,2 — Uj — Uj+i)) 


67 


3220734008 


max (min (0, 1 — Uj — m^+i) , min (mj-2 + Uj-i — 1,2 — Uj — Wj+i)) 


95 


3487613920 


max (min max (0, 1 — Uj — Uj+i)) , min {uj-2 + Uj-i — 1, 1 — Uj)) 


98 


3489185848 


max (min (0, 1 — Uj — Uj+i) , min {uj-2 + Uj-i — 1, 1 — Uj)) 



m 


N{f) 




2 


2881005752 


max (min {Fj,Fj-2 + 1) , Fj+2 - 1) 


13 


3099375756 


min {Fj+i, max {Fj,Fj-2 + l,Fj+2 - 1)) 


14 


3099506818 


min (Fj+2, max {Fj,Fj-2 + 1, Fj+2 - 1)) 


15 


3099572352 


min(Fj - l,^;_2,i^j+2 -2) 


53 


3202581216 


min (max (Fj_2, min {Fj,Fj+2 - 1)) , max (-Fj-2 + 1, -Fj+2 - 2)) 


63 


3219162080 


min (max (F^-i , min {Fj,Fj+2 - 1)) , max {Fj-2 + 1, -Fj+2 - 2)) 


67 


3220734008 


min (max {Fj,Fj+2 - 1) , max (F,-_2 + 1, -Fj+2 - 2)) 


95 


3487613920 


mill (max {Fj^i, iiiiii {F.j, Fj+2 — 1)) • niax {Fj^2 + 1- Fj+i — 1)) 


98 


3489185848 


miu (max (Fj,Fj+2 - 1) . max {Fj^2 + 1. -Fj+i - Ijj 
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m 


Alt f\ 




z 


0CC1 nnf;'7f;o 


'■ — 7 7 ^ \ 

min(niax(0, Xi—i + 2), Xi^i — 2) 




oUvvoi 0(00 


max(a;i — 1, min(0, Xi-i + 2), Xi+i — 2) 


1 A 


ouyyouooio 


niax(xi — 2, niiii(0, Xi—i -\- 2), i?^7-t-i — 2) 


15 


3099572352 


min(xi+i , + 2, .x.,+2 - 2) 


53 


3202581216 


max(min(a;i + 2, max(0, Xj+i — 2)), min(a;i_i + 2, a;i+2 — 2)) 


63 


3219162080 


max(min(xj + 1, max(0, Xi+i — 2)), min(a;i_i + 2, Xi+2 — 2)) 


67 


3220734008 


max(inin(xi, Xi+i — 2), min(xi_i + 2, Xi+2 — 2)) 


95 


3487613920 


max(inin(xi + 1, max(a;j, x^+i — 2)), niin(xj_i + 2,a;i+i — 1)) 


98 


3489185848 


max(min(a;j, Xi+i — 2), min(a;j_i + 2,Xi+i — 1)) 



Acknowledgments 

This work was supported by JSPS KAKENHI Grant Number 22560068. 
References 

[1] Wolfram S 1986 Theory and Applications of Cellular Automata (Singapore: World Scientific) 
[2] Wolfram S 2002 A New Kind of Science (Champaign: Wolfram Media) 
[3] Boccara N and Fuks H 1998 J. Phys. A: Math. Gen. 31 6007 

[4] Tokihiro T, Takahashi D, Matsukidaira J and Satsuma J 1996 Phys. Rev. Lett. 76 3247-50 

[5] Matsukidaha J, Satsuma J, Takahashi D, Tokihiro T and Torii M 1997 Phys. Lett. A 225 287-295 

[6] Nishinari K and Takahashi D 1998 J. Phys. A: Math. Gen. 31 5439-50 

[7] Takahashi D, Matsukidaira J, Hara H and Feng B 2011 J. Phys. A: Math. Gen. 44 135102 

[8] Hattori T and Takcsuc S 1991 Physica D 49 295-322 
[9] Matsukidaira J and Nishinari K 2003 Phys. Rev. Lett. 90 088701 
[10] Matsukidaira J and Nishinari K 2004 Int. J. Mod. Phys. C 15 507-15 



